Scaling Properties of Fidelity in Spin-one Anisotropic Model 
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By means of the density matrix renormalization group technique, the scaling relation of the fi- 
delity susceptibility proposed recently is verified for the spin-one XXZ spin chain with an on-site 
anisotropic term. Moreover, from the results of both the fidelity susceptibility and the entangle- 
ment entropy, the critical points and some of the corresponding critical exponents are determined 
through a proper finite-size scaling analysis, and these values agree with the findings in the litera- 
ture. Thus our work provides a numerical support of the use of the fidelity in detecting quantum 
phase transitions. 
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Quantum phase transitions (QPTs)j 1 - driven by purely 
quantum fluctuations, are characterized by the dramatic 
changes in the ground state of a many-body system as 
the controlling parameters in the system Hamiltonian are 
varied across critical points. Due to latest advances in 
quantum information science^ people attempt to char- 
acterize QPTs from the perspective of quantum informa- 
tion. One of the well-studied aspects is to explore the role 
of quantum entanglement in identifying QPTs.— In par- 
ticular, as a bipartite entanglement measure, the entan- 
glement entropy of a block of length I for one-dimensional 
systems is shown to exhibit qualitatively different scaling 
behaviors at and off criticality ^ 5 ' 6 ' 7 i 8 ' 9 ' 10 i u i 12 ^ The en- 
tanglement entropy saturates to a finite bound as the 
length I increases for noncritical (gapped) systems^^ 
whose value can vary for different boundary conditions. 7 
However, the entanglement entropy increases logarithmi- 
cally for critical (gapless) systemsj 4 ' 8 ' 9 ' 10 ' u i 12 i 13 By using 
conformal field theory, a universal scaling is expected at a 
quantum critical point, and its expression depends again 
on the boundary conditions. Thus the divergent charac- 
ter of the entanglement entropy in the finite-size scaling 
can faithfully indicate the existence of the critical points 
for one-dimensional systems. 

In the last few years, the ground-state fidelit y 14 : 15 (and 
its second derivatives, the so-called "fidelity susceptibil- 
ity"—) , another concept emerged from quantum informa- 
tion science, attracts much attention on their application 
to the analysis of QPTsi^Ol^dl 2 ^^^^ 2 ^^ 

As illustrated before in several concrete models, it seems 
that the singularity in the fidelity susceptibility can be an 
effective tool in detecting critical points. Quite recently 
general scaling analyses of the fidelity susceptibility are 
proposed i 23 : 24 As explicitly shown in Ref. [HI, the fidelity 
susceptibility S must be bounded above in the thermody- 
namical limit for noncritical (gapped) systems containing 
only local operators. However, for critical (gapless) sys- 
tems of finite size L, it fulfills scaling relations 
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where d is the spatial dimension, z is the dynamic expo- 
nent, and Ay is the scaling dimension of the transition- 
driving term in the Hamiltonian. This result implies that 



the QPTs at those critical points with Aq < can be 
detected by the power-law divergent behaviors in S. 

In this paper, the spin-one XXZ spin chain with a 
uniaxial single-ion anisotropic term is investigated, and 
we focus our attention on the verification of the predicted 
scaling behavior of the fidelity susceptibility in Eq. |T]). 
It is known that, while the QPTs can in principle be 
unveiled by the knowledge of the entanglement entropy 
and the ground state fidelity, they are usually difficult to 
be calculated due to the lack of knowledge of the exact 
ground state wavefunctions. Although numerical exact 
diagonalization can always be employed to evaluate the 
entanglement entropy and the fidelity for small systems, 
this method may not be able to demonstrate the scaling 
behaviors because of finite-size effects. Thus we need 
to resort to the density matrix renormalization group 
(DMRG) technique 2 ^ for the calculations for systems of 
large sizes. In the present work, both the entanglement 
entropy and the fidelity susceptibility are evaluated by 
means of the finite-system DMRG technique under open 
boundary conditions for system sizes up to L = 160. In 
our DMRG calculations, up to 300 states per block are 
kept and five DMRG sweeps are performed for the trun- 
cation error being about 10 -10 . We find that develop- 
ing peaks do appear in both measurements, which signal 
precursors of the QPTs. Applying a proper finite-size 
scaling analysis, the proposed scaling relation in Eq. ([1]) 
is confirmed numerically. Besides, the critical points in 
the thermodynamic limit and some of the correspond- 
ing critical exponents are determined through a proper 
finite-size scaling analysis, and these values agree with 
the results in the literature. Moreover, the results com- 
ing from both the entanglement entropy and the fidelity 
susceptibility are consistent each other. This implies that 
both measurements are equally suited for revealing QPTs 
and pinning down the critical points in the present case. 

The Hamiltonian for spin-one XXZ spin chains of L 
sites with an on-site anisotropic term is 

j=i i=i 

(2) 

where S" (a = x, y, z) are the spin-one operators at the 
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j-th lattice site. A and D parametrize the Ising-like and 
the uniaxial single-ion anisotropies, respectively. The 
full phase diagram consists of six different phase d 29 ! 30 
(see Refs. [3lll32ll33l for recent numerical determinations). 
On the A > half-plane, there consists of three phases 
known as the Haldane, the large- .D, and the Neel phases. 
All these three phases show a nonzero energy gap above 
the ground state. Between these phases, various types 
of phase transitions take place. There are a Gaussian 
transition between the Haldane and the large- 1? phases, 
and an Ising transition between the Neel and the Hal- 
dane phases. These two transitions merge at a tricrit- 
ical point A ~ 3.20 and D ~ 2.90^-^ where the Hal- 
dane phase disappears and the Neel-large-D transition 
becomes first order. Here we consider only the Gaussian 
and the Ising transitions at A = 1. In this case, it has 
been found that, as D is decreased from a large value, 
one first meets a Gaussian transition from the large- D to 
the Haldane phases at the critical point D c ~ 0.99, and 
then an Ising transition from the Haldane to the Neel 
phases at D c ~ -0.31j2!»2a 

For the convenience of the following discussions, some 
details of the Ising and the Gaussian transitions are re- 
viewedi Both transition lines are of second order with a 
dynamic exponent 2=1. Nevertheless, the former is de- 
scribed by a conformal field theory (CFT) with a central 
charge c = 1/2, while the latter by a c = 1 CFT.— More- 
over, their singular behaviors with the universality class 
of the transition, i.e., critical exponents, can be different. 
For the Ising transition, it is known that the correlation 
length critical exponent v = 1 and the scaling dimen- 
sion of the transition-driving term in the Hamiltonian 
Ay = l. 1 However, for the Gaussian transition between 
the Haldane and the large- 1? phase, it is found that the 
low-energy effective continuum theory can be described 
by the sine-Gordon mode l 32 ! 34 
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where n and $ are the conjugate bosonic phase fields, 
and a is a short-distance cut-off of the order of the lattice 
spacing. The coefficient /i oc (D — D c ) in the vicinity of 
the critical point D c for a given A, and thus becomes zero 
at the transition point. The value of the Luttinger liquid 
parameter K varies continuously between 1/2 and 2 along 
the critical line. We note that all the scaling dimensions 
and the critical exponents are determined by a single 
parameter K. Consequently, they change continuously 
along the critical line. From the sine-Gordon theory^ 
it is found that the critical exponent of the correlation 
length v = 1/(2 — K) and the scaling dimension Ay = K 
for the transition-driving term cos(v47rX ( I>). 

In the following, our DMRG results are presented in 
orderi 3 ^ The findings of the fidelity susceptibility S(D) 
and the ground state fidelity T(D, D + 8) are shown in 
Fig. ffl The fidelity susceptibility, or the second derivative 




FIG. 1: (Color online) Fidelity susceptibility S(D) for the 
spin-1 XX Z spin chain in Eq. (J2Jl as functions of D for various 
sizes L with A = 1. Inset shows the fidelity T{D, D + 8) 
as functions of D for the corresponding sizes. Here we take 
5=1(T 3 . 



of the fidelity, is calculated b y 18 ' 20 

2[l-F{D,D + 8)} 



S(D) = lim 

5^0 
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where the ground-state fidelity (or the modulus of the 
overlap) is given b y 15 i 37 



F(D,D + S) = \{9 (D)\9 (D + S))\ 



(5) 



with \^ (D)) and |^'o(-C ) + 8)) being two normalized 
ground states corresponding to neighboring Hamiltonian 
parameters. In our calculations, 8 = 10~ 3 is used. As 
shown in the inset of Fig. [T] drops in the ground state 
fidelity are observed, which signal precursors of the Gaus- 
sian and the Ising transitions in the model under consid- 
eration. The drops in !F(D, D + 8) at the right-hand side 
show the Gaussian transition, while those at the left-hand 
side give the Ising one. Further evidences for indicating 
QPTs are provided by the results of S(D). As seen from 
Fig. HI the maximum values <S max in the fidelity suscepti- 
bility grow with increasing size, and thus indicate diver- 
gence in the L — > oo limit (see also Fig. @] below). From 
the scaling analysis in Ref. l23j, these divergent behaviors 
in iS must imply the appearance of the QPTs. Applying 
the finite-size scaling, the critical points D c in the ther- 
modynamic limit can be determined from the locations 
D max (L) of the local maxima in S(D) on a size-L system 
(see Fig. |3] below). 

As mentioned before, the divergent character of the 
entanglement entropy can also show the existence of the 
QPTs. Thus the entanglement entropy is evaluated for 
comparison. Here we consider the entanglement entropy, 
or the von Neumann entropy of the reduced density ma- 
trix pn.(D) of the right-hand block of L/2 contiguous 
spins 



£(D)=-Tr{p R (D)log 2 p R (D)] 



(6) 
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FIG. 2: (Color online) Entanglement entropy £ (D) for the 
spin-1 XX Z spin chain in Eq. (J2J) as functions of D for various 
sizes L with A = 1. 



Our DMRG results are shown in Fig. [5J It is found that, 
far away from the critical points, £ (D) has no size depen- 
dence, as expected for the gapped phases. Nevertheless, 
two peaks develop as size L increases. Again, these peaks 
indicate the existence of the Gaussian and the Ising tran- 
sitions, and the corresponding critical points D c can be 
deduced from the locations D max {L) of the local maxima 
in £(D) on a size-L system, as discussed below. 

According to the finite-size scaling theory,— one has 



\D max (L) - D c \ cxL 



-l/u 



(7) 



where D c is the critical point in the thermodynamic limit 
and v is the critical exponent of the correlation length. 
Thus D c can be determined by an extrapolation pro- 
cedure. The results for the Haldane-Large-Z? and the 
Haldane-Neel transitions are shown in Fig. [3J For the 
Gaussian transition between the Haldane and the large- 
ly phases, both extrapolations give D c ~ 0.97 as shown 
in the top panel of Fig. [3] The critical exponent of the 
correlation v ~ 1.42 for the data of Z) max (L) obtained 
from £, while v ~ 1.45 for those from S. Because of the 
relation v = 1/(2 — K), the Luttinger liquid parameter 
K = 1.30 [K = 1.31) for the data related to £ (S). We 
find that the values obtained from the measurements of £ 
and S agrees each other, and they are consistent with the 
previous findingS ) 32 i 39 where D c ~ 0.99 and K ~ 1.328. 
For the Ising transition between the Neel and the Haldane 
phases, D c ~ —0.31 for both extrapolations as shown in 
the bottom panel of Fig. [3] The critical exponent of the 
correlation v ~ 0.90 for the data of Z? max obtained from 
£, while v ~ 1.05 for those from S. Again, the value 
of D c agrees with the previous result (D c ~ — 0.31) j 33 i 39 
and our findings of v are consistent with the theoreti- 
cal prediction (y = 1) for the Ising transition. From the 
above discussions, we find that both the entanglement 
entropy and the fidelity susceptibility are equally suited 
for revealing the critical behaviors in the present case. 
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FIG. 3: Finite-size scaling of D m ax versus L~ x l v . The full 
lines are least square straight line fits for sizes with L > 100. 
Top: the Haldane-Large-D transition, where v ~ 1.42 (y ~ 
1.45) for those D ma x's corresponding to the local maxima in 
the curves of £ (S). Bottom: the Haldane-Neel transition, 
where v ~ 0.90 (y ~ 1.05) for those D ma x's corresponding to 
the local maxima in the curves of £ (S). 



To verify the predicted critical scaling behavior of the 
fidelity susceptibility in Eq. fT}, the values S ma , x (L) of the 
local maxima for various sizes L are plotted in Fig. [4j It 
is found that our data do fulfill the scaling relation in 
Eq. ((I]), where A Q = -0.33 (i.e., Ay = 1.34) for the 
Gaussian transition and Aq = —0.89 (i.e., Ay = 1.06) 
for the Ising one (d = 1 and z = 1 are assumed here). 
The value of Ay for the Ising transition agrees with the 
predicted one, Ay = 1. Since Ay = K for the Gaussian 
transition, the Luttinger liquid parameter K determined 
by the present finite-size scaling agrees with the previous 
findingsSS, and those determined by the critical exponent 
v coming from the scaling in Fig. [3] Thus the fact that 
a single parameter K controls all the critical exponents 
for the Gaussian transition is confirmed by our numerical 
results. 

In summary, the general scaling analysis of the fi- 
delity susceptibility proposed in Ref. [23| is verified by the 
present DMRG calculations for the model of Eq. The 
critical points of the Gaussian and the Ising transitions, 
as well as some of their critical exponents, are determined 
from the perspective of quantum information. We note 
that, as seen from Figs. [3] and [H data for systems of 
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FIG. 4: The log-log plot of 5 max for various sizes L. The full 
lines are least square straight line fits for sizes with L > 100. 



smaller sizes can deviate from the fitting lines obtained 
from the data for those of larger sizes (say, L > 100). 
Therefore, to avoid the finite-size effects and to unveil 
the correct scaling behaviors at the critical points, cal- 
culations for systems of large enough sizes are necessary. 
From our DMRG calculation for systems of large sizes, 
we conclude that the fidelity susceptibility and the en- 
tanglement entropy can have similar predictive power for 
revealing QPTs. 
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We note that the DMRG calculations in Refs. and H3 
are under the periodic boundary conditions, rather than 
the open boundary conditions used in the present work. 
The small discrepancy of our results from theirs may be 
due to the different boundary conditions employed. 



